Half-knot in the spinor condensates 
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We present an exact solution to the stationary coupled nonlinear Gross-Pitaevskii equations which 
govern the motion of the spinor Bose-Einstein condensates. The solitonic solution is a twisted half- 
skyrmion in the three-dimension (3D) space. By making a map from and the Cartesian coordinates 
to the toroidal coordinates, we demonstrate it is a linked half-unknot with a fractional Hopf charge. 
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Topological objects appear in a variety of branches of 
physics from the liquid-crystal colloids [T], the classical 
field theory[51 [3], and the nonlinear optics[31 [S], to the 
condensed matter physics [6l [7] and the modern universe 
theory^S . In the dilute atomic gas, quantized vortices 
and skyrmions have been widely studied and experimen- 
tally observed. In particular, the spinor Bose-Einstein 
condensation (BEC) whose spin direction can change dy- 
namically has provided an ideal pilot to model and ex- 
plore the rich family of topological excitations [9} llOj. 

The Hamiltonian of the spinor BEC has the SO{3) x 
U(l) symmetry. In contrast to '^He superfluid where 
only the gauge symmetry is broken, the dilute atomic 
spinor BEC can also break the spin-rotational symme- 
try. A large amount of theoretical works has been carried 
out to study the nontrivial topological structures such 
as the fractional vortices, the skyrmions, the monopoles 
and the knots as well [TT| - [T7] . It is known that two 
kinds of spin textures, the Anderson- Thoulouse skyrmion 
and Mermin-Ho skyrmion, exist in 2D multi-component 
BECs. In 3D spinor BECs, monopoles, 3D skyrmions, as 
well as knots are theoretically predicted and numerically 
simulated. 

According to the homotopy theory, the topological ex- 
citations are classified by the broken symmetry of the 
corresponding ground states [9l [TOl [18]. They are char- 
acterized by the topological invariants. In the antiferro- 
magnetic spinor BEC, the order parameter manifold is 
M = X U{l)/Z2. Here U{1) denotes the manifold 
of the superfluid phase (f> and S"^ is a two-dimensional 
sphere. The homotopy groups yield 7r„(M) = Z {n = 
1,2,3). Singular line defects (vortices) and point defects 
(monopoles) are determined by the first and the second 
homotopy classes, 7ri(M) and 7r2(M), respectively [13, 14]. 
The monopoles are sometimes called singular Skyrmions. 
The nontrivial third homotopy class tt^ (M) implies topo- 
logical objects of knots. The corresponding topological 
invariant is called the Hopf charge which counts the num- 
ber of times the 2D sphere covered by the 3D sphere. 

The mean-field order parameter of the spinor BEC 
is governed by the coupled Gross-Pitaevskii equations 



(GPE). We present the exact solitonic solutions to the 
stationary GPEs for the F — 1 BEC in a uniform mag- 
netic field. To our knowledge, this is the first analytical 
result to the coupled nonlinear equations in the 3D sys- 
tems. Previous works are mainly in lower-dimensional 
space or/and by means of numerical simulations, in which 
the rich possibility of intriguing topological structures 
are often elusive[3S]. Attempts for seeking analytical so- 
lutions to the coupled GPEs are usually frustrated by 
two obstacles: one is the nonlinear density-density in- 
teractions; the other is the nonlinear spin-spin interac- 
tions. We propose a method to decouple the equations 
by reassembling the spatial dimensions and the hyperfine 
states in the nonlinear density-density interactions. The 
complexity of the spin-spin couplings are surpassed by de- 
manding that the hyperfine states satisfy ijj-i = — V'l and 
■0Q = ipo, which yield the average of the spin-polarization 
|F| = 0. This vanishing spin-polarization state is most 
easily realizable in the antiferromagnetic BECs. The soli- 
tonic solution represents a twisted half-skyrmion in the 
three-dimension (3D) space and is demonstrated to be a 
linked half-unknot with a fractional Hopf charge. 

The mean-field Hamiltonian reads, 

dr J2 CJ-^V2-pm + (7m2 + y(r)]V'm 

m=0,±l 

+ ^ J dr[cop^ + C2F% (1) 

where the terms associated to p and q represent the linear 
and the quadratic Zeeman effects, respectively. p(r) — 
J2m. lV'm(r)P is the total density of the condensate. The 
spin-polarization vector F — J2m ''Pm^mnipn/ P, with f the 
spin matrices. The nonlinear coupling constants cq — 
((70+252)73 and C2 = (52— 5o)/3, where relating to the 
s-wave scattering length of the total spin-F channel as 
gp = ^Trh^ap /M . By choosing the units of ?i = M = 1, 
the stationary GPEs in a uniform external field {V{y) = 
0) are written as 
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[-2 V -P + q+ (co + C2)(|V'i| + IV^ol") + (co - C2)|7/^_ir]Vi + C2Vod, 



MV'O 



+(co + C2)(|^iP + l^-iP) + colV'oPlV'o + 2c2V'SV'iV'- 



V^+P + q+ (co + C2)(|^_l|' + IVoH + (co - C2)|7/'-l|^]V'-l + C2V'oV'l* 



(2) 



r 



The 3D solitonic form of solutions is constructed as 



— ^[tanh(A;x)sech(fc?/) — isech(fca;)tanh(fcy)]e 



sech(fcx)sech(fci/) 
[tanh(/ca;)sech(A;?;) + jsech(A:x)tanh(fcy)]e* 



_(3) 

where the periodic boundary condition in the z-direction 
requires kz = 2mr/Lz with n the number of modes. For 
simphcity, we have chosen = ky = k. The total den- 
sity of the condensate p{r) = sech^(fca;) + sech^ (ky) is 
uniform in the z-direction, which endues the Eqs.([2]) be 
decoupled and the variables be separable according to the 
dimensions. This solution has solitonic form in the x — y 
plane. The parameter relations are straightforwardly ob- 
tained by substituting the state ([3| into the Eqs.([2]) as 
Co " 



the quadratic Zeeman energy plays the role of balancing 
the chemical difference between the hyperfine states. 

The 3D order parameter can be rewritten as ^(r) — 
\J /3(r)x(r), where the x(i") is a normalized spinor, 
X^(r)x(r) = 1. Instead of the vanishing spin- 
polarization, we consider a real unit vector d ~ 
{dx,dy,dz)'^ , by which the spinor x{^) is represented 
aslMESI 



x(i-) = 



/ -d^ +idy \ 

dz 



-fc^, p = 0, g = — and /i — —k^. We note that The d(r) vector is then given by 



(4) 



d(r) 



(tanh(fca;)sech(fci/) cos(A;2 2:) — sech(A;x)tanh(A;2/) sm{kzz) 
tanh(fcx)sech(A;y) sin(fczz) -I- sech(fca;)tanh(fcy) cos{kzz) 
sech(fca;)sech(fcj/) 



(5) 



r 



In quantum magnetism, d is also called the Neel vector 
which corresponds to the staggered magnetization of the 
antiferromagnetic state. Since |d| = 1 and d^ > every- 
where, it forms the surface of a semi-sphere in the order 
parameter manifold which is denoted as S'J 
subscript 1/2 indicates half of a sphere 



^^2- Here the 



To reveal the topological structure of this state, we 
first examine the 2D situation which is fulfilled by sim- 
ply setting z = in the state ([s]). Figure 1(a) and (b) 
display the density distributions of the hyperfine states 
tpi (or and tpo, respectively. It shows that the hy- 

perfine state ipi {ip-i) forms a line-shape soliton, whereas 
■00 forms a 2D point-shape soliton. The phase distribu- 
tion of the ipi (V'-i) in Fig. 1(c) reveals a vortex (anti- 
vortex) at the center, where the tpo component resides in. 
Therefore, the spinor BEC forms a coreless solitonic vor- 
tex. Fig. 1(d) displays the spherical harmonic representa- 



tion of the spinor condensate density which is expressed 
by |^'(r,6l„(/),)|2 = |Em=o,±iV'm(r)Yi™(6's,(/>s)P, where 
{6s,<j)s) are the coordinates in the spin-space fgl. It ex- 
hibits a semi-hedgehog structure. We demonstrate this 
structure is a half-skyrmion as follows. 

In 2D geometry, the vector d(r) — 
—J= (tanh(fca;)sech(fc2/) , sech(A:x)tanh(A:?/) , sech(fca;)sech(fcy))'^ 

still represents the surface of a semi-sphere. The vecto- 
rial distribution of the d field is shown in Fig. 2(a). The 
vector points upward at the center and gradually twists 
in the 2D plane at large distance. This texture represents 
a half-skyrmion which covers exactly half a sphere [5U]. 
It is also called a Mermin-Ho Skyrmion or merlon which 
has been studied in the condensed matter physics such 
as the quantum Hall effects [12] and high temperature 
superconductivity [5D]. It has also been created in the 
superfluid ^He-A in a rotating cylinder [55] and in 
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FIG. 1: (Color online) The 2D density distribution of the state 
(|3| for (a) Tpi (or ip-i) a-nd (b) tpo- (c) The phase distribution 
for ^-1. (d) The spherical harmonics representation of the 
spinor density. 




FIG. 2: (Color online) (a) The 2D texture of the d field, (b) 
The topological charge density g(r). The total topological 
charge is Q = 1/2. 

dilute atomic BEC by adiabatic deformation of the 
magnetic trap[TT]. Fig. 2(b) shows the topological charge 
density q{r) = ^d - {dxd x dyd), which characterizes the 
geometric distribution of the spin texture [21]. The total 
charge or the Pontryagian index Qp = / g(r)dr = +1/2 
is a topological invariant |231 [21] . This is the reason 
why it is called a half-skyrmion. The anti-half-skyrmion 
solution with Qp = —1/2 can be obtained by simply 
taking the complex conjugacy of the state ([s]). 

Now we further explore the 3D topological structure 
of the state The periodic boundary condition in the 
z-direction implies that it essentially forms a twisted vor- 
tex ring. Figure 3 display the vectorial distribution of the 
d(r) vector in the 3D space. The vector points from up- 




FIG. 3: (Color online) The 3D texture of the d field. 

ward at z = to inplanc at large distance while twists 
around the z-direction. The state that has a vanishing 
average spin-polarization implies it is locally degenerate 
[14]. The order parameter is invariant under an arbi- 
trary rotation about d, i.e., exp[— id • ft/i]^' — "if, where 
(/) is the rotation angle. The manifold of the unit d vec- 
tor in the 3D space is similar to that in the 2D space. 
Since x(r) is independent of any particular direction, it 
represents a mapping from a compactified 3D domain 
(r) into a 2D semi-sphere in the order parameter space 
^ 3^2- The preimage of a point on the tar- 

get corresponds to a closed loop in S^. It can be 
characterized by a topological invariant called the Hopf 
charge, in integral formula which counts the number of 
times the two-sphere is covered by the three-sphere 




FIG. 4: (Color online) Schematic description of the mapping 
Q from a cylinder to a torus. The symmetric axis of the 
cylinder is mapped as the core of the torus. 
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In order to get a more intuitive view of the topolog- 
ical structure, we make a mapping from the Cartesian 
coordinates to the toroidal coordinates by 

x\ / i\ / {x + R)cos{z/R)\ 
y j ^ iy j = {{x + R) sin(z/i?) j , (6) 

where R — Lz/2tt defines the size of the torus. This 
mapping transforms a cylinder with periodic boundary 
condition into a torus, as schematically shown in Fig. 4. 
|r| — R indicates the position of the toroidal core. In 
terms of vector d it twists n-times around the core before 
joining the ends, implying it is a twisted vortex ring. Vor- 
tex rings have been studied experimentally in superfluid 
liquid helium [8, 9] and in Bose- Einstein condensates [10]. 
Numerical simulations have revealed that superfluid tur- 
bulence contains linked vortex lines [11]. Since we have 
the surface of a semi-sphere instead of a full sphere in 
the order parameter manifold, there are two types of dis- 
tinct vortex rings: one is a twisted tube-ring; the other 
is a twisted Mobius strip, depending on the choice of the 
points on the S'^^^ ■ Figure 5 schematically shows the two 
distinct mappings for the twisted tube {d^ = 0.98) and 
the twisted Mobius strip {dx = —0.98). The color on the 
tube represents for the value of arctan(c?j^/dj;) while on 
the Mobius strip the value of arctan((ij^/(i2). It helps to 
recognize the twisting feature of the tube or the strip. 
We mention that the color on the semi-sphere is just 
schematic and has no direct relevance to the color on 
the tube or the strip. 




FIG. 5: (Color online) A loop (string) on the Sf/2 corresponds 
to a twisted tube (Mobius strip) in the real space. The color 
indicates the twisting number. 

In knot theory this twisted vortex ring defines an topo- 
logical unknot, the simplest knot-like structure [Hinj. The 




FIG. 6: (Color online) Two loops or two strings in the order 
parameter manifold correspond to two linked tubes (a) or 
two linked Mobius strips (b) in the real space. The twisting 
number equals to the linking number of the tubes or strips, 
e.g., (c) and (d) One tube and two Mobius strips for n — 1 
and n = 2, respectively. 



knot can be characterized by the Hopf charge which is 
defined by [7] 

Qn = ^ J d.^^£i.jkJ^ij^k, (7) 

where J-ij — diAj — djAi — d • [did x djd) is the strength 
of the gauge field and the vector A(r) defines a connec- 
tion in the order parameter space. If the d field has Hopf 
charge Qh, the tube or strip twists Qh times. Further- 
more, two loops corresponding to the preimages of any 
two distinct points on the target S'^^^ ^il^ be linked Qh 
times. It is directly calculated that the state ^ has 
a fractional Hopf charge of Qh ~ where the factor 
1 /2 comes from the semi-sphere instead of the full-sphere. 
Hence we call it the half-knot so as to differentiate it from 
the usual integral number of Hopf charge. 

The topological structures are shown in Fig. 6. Fig. 6(a) 
is two linked twisted tubes which corresponds to two 
closed contours on the Sl^^ while Fig. 6(b) is two linked 
Mobius strips which corresponds to two open strings for 
dx = 0.98 and dy = 0.98, respectively. Fig. 6(c) and (d) 
are examples for three linked loops for n — 1 and n — 2, 
respectively. The linking number is equal to the twisting 
number of the tubes or strips. 

In summary, we have constructed an exact solution 
to the coupled nonlinear GPEs in the 3D space. We 
revealed that it represents a solitonic knot with a half- 
integer topological charge. Our method may provide a 
clue to analytically solve other coupled nonlinear equa- 
tions. Some relevant works will be published elsewhere. 
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